Photoemission in the system of linear chains : 
apphcation to PrBa2Cu307 and La2 x yNdySrxCu04 
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Photoemission in the system of linear charged chains has been studied with the dipole matrix 
incorporated, from which its dependencies on the photoelectron momentum (k), photon polarization 
(E), and photon energy (u) can be explored. The used model is the three-dimensional array of non- 
interacting chains, which is so simple as to allow an analytic approach. Motivations of the study are 
for the doped CuOs chain in PrBa2Cu307 and the doped static stripe phase in La2-x-yNdySrxCu04. 
The one-dimensional dispersion exhibiting spin-charge separation and its dependence on the mo- 
mentum perpendicular to the chain and its E-dependency as well are discussed in PrBa2Cu307. 
For La2-x-yNdySrxCu04, the anomalous spectral distribution formed by two sets of stripes per- 
pendicular to each other is investigated. The geometric effects led to by the dipole transition matrix 
including the interference of photocurrents from different chains are found to change the simple 
one-dimensional feature drastically. We find these changes are consistent with experiment for the 
chain system PrBa2Cu307, but less satisfactory for the stripe phase in La2-x-yNdySrxCu04. This 
means that in the stripe phase much two-dimensional characters still exist unlike the chain system. 
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I. INTRODUCTION 



Angle-resolved photoemission spectroscopy (ARPES) 
has become the most significant experimental tool for a 
search of the electronic structure, especially in the high 
temperature superconducting cuprates p| . In many the- 
oretical and experimental interpretations, the ARPES 
data is often compared to the one-electron spectral func- 
tion A(k, Lo) or vice versa. For a more satisfactory de- 
scription of the photoemission, i.e., if one is interested 
in various parameter dependencies or the line shapes of 
quasi-particle peaks, one must go beyond the spectral 
function . For it, we may need to consider the photoex- 
citation process governed by the dipole transition matrix 
|^,|| and the extrinsic losses of the photoelectron on its 
way out to and through the surface |^,^. In this paper, 
we give the study of transition matrix effects in the pho- 
toemission, while the extrinsic losses from the dynamical 
scattering are not our concern because their contributions 
are in many cases smaller by orders of magnitude. 

The target systems are composed of charged linear 
chains, i.e., one-dimensional metals. One-dimensional 
metal has been known for the extreme realization of 
electron correlation effects in that there are an exotic 
spin-charge separation and only collective modes, not 
the single-particle excitation [Q. It has opened an- 
other paradigm different from Fermi liquid and is called 
Tomonaga-Luttinger liquid. In an actual situation, a 
one-dimensional metallic state is very unstable toward 
an insulating state or very marginal between metallic 
and insulating states through charge and spin ordering. 
Anyhow these one-dimensional chain systems have at- 
tracted considerable attentions especially since the dis- 



covery of high Tc superconductors. Among the high 
Tc cuprates, YBa2Cu307 (YBCO) and its family com- 
pounds have CuOs chain structure as well as Cu02 
plane. PrBaaCuaOy (PBCO), among the YBCO fam- 
ily cuprates, does not support superconductivity un- 
like other cuprates, which is suggested due to the hole 
depletion in Cu02 plane IqM- Additional holes are 
thought to be doped in the CuOs chain, not in the 
Cu02 plane. In such a sense, PBCO can probably be an 
archetype of Tomonaga-Luttingcr-type one-dimensional 
system, which has revived much interests in the study 
of metallic chain system within the study of high Tq 
cuprates (PBCO actually has a band gap opening possi- 
bly due to charge ordering in the chain IQ]). 

On the other hand, another feasible usefulness of the 
metallic chain system in high Tc superconductor physics 
is in the so called stripe phase observed in Nd-substitutcd 
La2_x-yNdySrxCu04 (Nd-LSCO) |ll|,|li- The stripe 
phase is due to the formation of an ordered array of 
charged stripes which are also antiphase domain walls 
between antiferromagnetic ordered spins in the Cu02 
planes, which gives the periodicity-doubling for spins 
compared to charges. 

Recent ARPES experiments for PBCO ||l^ have shown 
actual evidence for one dimensionality immersed in, 
but distinguishable from contributions from the two- 
dimensional Cu02 planes. In experiments, two kinds 
of dispersive band features are observed; one is from 
the undoped Cu02 plane and the other is the doped 
CuOa chain. Especially in the one-dimensional-like fea- 
tures, the separation between charge and spin excitation 
is also suggested. It is a stimulating result of Tomonaga- 
Luttinger liquid that the low energy physics is dominated 
by uncoupled collective modes of charge and spin exci- 



tations, also called holons and spinons. In experiments, 
it is also found that, in addition to the spin-charge sep- 
aration, one-dimensional characters (particularly spinon 
parts) are strongly dependent on the photon polarization 
and at the same time on the momentum perpendicular to 
the chain. The results for PBCO should be also relevant 
to other one-dimensional insulator, SrCu02, which has 
a weakly coupled double Cu-0 chain. It is very intrigu- 
ing to find that the similar polarization and perpendic- 
ular momentum dependencies in ARPES are observed 
in SrCu02 as well [|l3|. These observations could be 
hardly explained by a single Tomonaga-Luttinger chain 
only, the additional effect probably needs incorporating, 
which motivated us to explore the dipole transition ma- 
trix related to the geometries of a chain system. 

The static stripe phases in Nd-LSCO also lead to an 
interesting puzzle on the electronic structures through 
ARPES observations ||lj]. Prior to that, the neutron 
scattering experiments compcllingly proposed that there 
should be the static stripe formation in Nd-LSCO and 
besides, the charged stripes be aligned with rotated by 
7r/2 between two adjacent planes ||l2|. The low bind- 
ing energy electronic structure actually looks more or 
less one-dimensional in the sense that the spectral weight 
distribution is roughly understood as a superposition of 
two perpendicular chains. However, the suppression of 
d-wave nodal signal along the Brillouin zone diagonal di- 
rection, particularly around the F point, could not be 
satisfactorily explained without additional physics. We 
suggest one necessary ingredient could be from the dipole 
transition matrix including the interferences from differ- 
ent stripes. Another noticeable point is that there have 
been no decisive evidences for Tomonaga-Luttinger liquid 
in the stripes phase unlike chains in PBCO. Therefore, 
in our study for Nd-LSCO, although we consider the ma- 
trix effects from one-dimensional system, we assume the 
simple Fermi liquid for a single stripe having well-defined 
Fermi point. Further, whether each stripe is a Fermi liq- 
uid or a Tomonaga-Luttinger liquid does not affect the 
conclusion of the present work because the calculated 
polarization and photon energy dependencies of ARPES 
spectra come from the matrix element effects reflecting 
the interference of photo electrons from different stripes. 

In this work, we account for the dipole matrix effects 
for the semi-infinite (z > 0) three-dimensional array of 
metallic chains, whose Green's function is considered. 
The slight photoelectron damping is taken into account in 
the simplest way. It is assumed that all the excitations 
including both charge and spin occur along the chain- 
direction, that is, the energy scale for the perpendicular 
excitation is very high. This means that we neglect the 
interchain (interstripe) interactions such as hybridization 
or Coulomb interaction, which may exist in the actual 
system and also the fluctuation of chain or stripe perpen- 
dicular to it, which could be important in a dynamical 
stripe phase. We think of two configurations. First, we 
simply model the three-dimensional array of chains for 
PBCO. On the other hand, for Nd-LSCO, we slightly ex- 



tend the model and consider the three-dimensional chain 
array perpendicularly crossed in adjacent planes. For 
both configurations, the model is analytically tractable. 
The resulting spectra is determined by complicated inter- 
ferences of chain contributions and thus the dipole ma- 
trix effects are geometric. The effects naturally give k- 
and E-dependent photoemission spectra, which is found 
consistent with the recent experiments. 

The outline of this paper is as follows. First we discuss 
the basic theory and the suitable starting point to ac- 
count for the dipole transition matrix in Sec. II. In Sec. 

III, we investigate the model for PBCO and discuss its re- 
sults. We can see how the spectral function of Luttinger 
liquid changes by the transition matrix. Such changes 
are also found quite consistent with experiments. In Sec. 

IV, we extend and constitute the model for Nd-LSCO 
and obtain the spectral weight distribution through the 
similar type of calculation. The calculation is also more 
or less consistent with the low energy weight distribution 
in experiments (here we assume a priori metallic stripes) 
in a limited case, which is from the interference among 
chains. In the final Sec. V, we present the concluding 
remarks and give the outlook. 



II. BASIC FORMALISM 

To account for the dipole transition matrix in the pho- 
toemission spectral intensity, we follow Pendrey's for- 
mula ||l5[ . Its benefit is to treat only one-particle Green's 
function without solving for any static stationary states. 
From Pendrey's the photoemission intensity /(k, E, uj) is 

I{k,E,uj) = --ImM{k,E,uj), 

n 

Af (k, E, oj) = mC+iE + uj)AG+{E)A^G-{E + w)|k), 

(1) 

where G* is the retarded or advanced one-particle 
Green's function and |k) the photoelectron state. The 
photoelectron state could be represented as a plane wave 
in the lowest approximation, even if the real photoelec- 
tron wave function may be damped through the scatter- 
ing with possible excitations in the solid, i.e., could be 
quite different from a plane wave except for in the very 
high energy region. In our study, assuming the solid oc- 
cupies z > space, we take into account the damping in 
the simplest way by taking the photoelectron wave as a 
plane wave in the in-plane direction and a damped wave 
in the z-direction. This would not be a very good rep- 
resentation due to the anisotropy inherent in the system 
even in the in-plane direction |1^ . We would expect and 
show, nevertheless, this could catch the essential physics 
as it should have. Then in Eq.(Q) we replace |k) by |k), 
which is given as 
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kz is a complex number, k^ = \/k1 + iV, and thus the 
photoelectron gets damped. It is worth noting that a 
damped photoelectron state is the formal solution of the 
one-electron Schrodonger equation 



(ek-Wo-S)lk) =0, 



(3) 



where E is the self-energy operator involving the 
photoelectron-solid interaction. Equation (g) can be un- 
derstood from the first approximation for its self-energy, 
Y, = —iTd{z). If we insert the plane wave states in Eq.(|]), 
the photoemission matrix M{k,E,uj) is dissolved as 

M(k,£;,cj)= /"dk'dk"dk"'dk""(k|G+|k')(k'|A|k") (4) 

x(k"|G+|k"')(k"'|At|k"")(k""|G^|k), 

where Gi and G2 denote the one-particle Green's func- 
tions at E and E + ui, respectively, and A ~ E • p, we 
would get, assuming the linear polarized photon, 

(k'|A|k")-(E.k')<5(k'-k"), k=(K,kz). (5) 



III. PHOTOEMISSION IN THE SYSTEM OF 
CHAINS : FOR PrBaaCuaOr 

The low energy features in ARPES for PBCO are ex- 
pected to be governed by the CuOa chain. In a for- 
mula unit, i?Ba2Cu307 (i?=rare earth) has pairs of Cu02 
planes between CuOs chain layers. Among those com- 
pounds, for the case of PrBa2Cu307, various experimen- 
tal [|| and theoretical studies j|] suggests that the Cu02 
plane should not be doped, but the CUO3 chain doped. In 
this study, we do not consider the undoped Cu02 plane 
since its contributions will be pushed to high binding en- 
ergies. 

The chains are along the y-direction (5-axis) . Then the 
formal Green's function of the system in the real space 
representation can be written as 



G±(r, 



j:j: 



E-e,±iO+ ' 



(6) 



where -07"" {^) ^re eigenstates of the system, '0™™ (r) = 
'4's{y)x{x: — C||m)x(z — c_Lm'). Here ^l^siv) is assumed to 
include both charge and spin degrees of freedom, while 
in the direction perpendicular to the chain the minimum 
excitation energy is assumed to be so high that x(x) or 
x(^) is taken the lowest localized one-dimensional state. 
Besides, we assume the overlap between the chains is very 
small. Then, if we rewrite Eq.(p|) in a more explicit way, 

00 
G±(r,r')=G±(2/-y') ^ x{x - c^\m)x* {x' - c\\m) 



X\z~cx_m - —)X \z - cj_m ~ ^), 



m'=0 



where we have assumed that z = is the surface and the 
first chain layer is positioned at z = ^. G^{y — y') is an 
ordinary single-particle Green's function in the Luttinger 
liquid including both charge and spin degrees of freedom. 
In the Luttinger model, the single-particle Green's func- 
tion and its property are already well known |17]| . This 
trivial factorization of the Luttinger Green's function is 
thanks to the complete suppression of the perpendicular 
excitation by the starting assumption. It should be also 
noted that, for the interchain distances, experiments give 
Ca^ Cb^ C|| and Cj_ ^ 3c|| ^. 

Now we see (k|G|k') is 



(k|G|k') = G{ky)S{ky - k'y) £ r{K)[r{k'z)Y 



^x{K)x{K) 



l5]<5(fc.-fci + G,| 



Gil 



, (8) 



where Gii = 2n7r/cj[ is the in-plane reciprocal lattice vec- 
tor due to the periodic lattice. For the explicit calcula- 
tion, we propose a suitable exponential form for x(a:) (or 



x(a;) = ^e-N/^ x(fc.) = 



TT 1 + kW '' 



(9) 



and additionally, 




l/a + ikz l/a + ikz_ 



X"(fc.)=x(fc.)e^'=^^'""+'"/'^. 



^{c±rn+c±/2) 



(10) 



Here we basically assume cy, c_l » a. In Eq.(p^), x"^{kz) 

is reduced into x™(fcz) when replacing kz by kz. Using 
Eqs.(||),(||), and ([lO|), we get the photoemission matrix 
M(]i,E,uj), aside from irrelevant constants, 



M{k,E,LU) = \G2{ky)\^G+{ky)\x{k,)\^ 



(11) 



m rn' rn" 



(7) 



/'dfc;dfc;'^^|x(fc. + G||)nx(fc. + Gj|)p 

G„ g;, 

x[exikx + G||) +eyky + ezk'^] 
X [e4kx + G\)+ eyky + Szk"] 

x{x"'{k'z)rx"''{K){x"''iK)rx"'"iK), 



where E = {ex^SyL^) and m, m! , and m" run from 
to +00. Equation (|ll| ) involves two types of summations 
with respect to G\\ . The summation for the reciprocal 
lattice vector can be done by use of 



M(k,i?,c.) - \G2{ky)\^Gt{ky)\x{K)\'\x{k.)?e 



--Im.fc^cj_ (■m+m ) ikzC± (m-—m ) 



m m' m" 



E 



e^'^ii^^ = ci 



00 

E 

rn— — 00 



(5(a; — CiiTTT,). 



(12) 



First, using Eq.(|l2|), let us note and define 

Gil 



TT ^[l + (fc,+G||)2a2]2 
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C|| 



'^e''=="=ii™[cji|m|/a+l]e-=iil"l/" 



X [£^C,{kx) £,{ra — m')S,{ra ~ m,") 
+£lklakxfi{m - m')^{m' - m") 

+el^{kxfCim - m')CK - m") 
+2exeykyS,{ka;)C{kx)^{m - m')l{m' - to" 
+eyezfcyC(fcx)^'f(TO - m!)C,{m' - to") 
+£yezfcy?(fcx)^C(w - m')^{m' - to") 

+£zea;C(^a;)C(fca;)C(™ " m')C,{m' - to") 

+£.e.C(^r.)C(fcr.)C(m - m')f(TO' - to")]. 



e(fc.), 



(13) 



and in the same way, 
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(fc,+G||)|x(fcx + G||)|" = — 



kx + Gil 



TT ^[l + (fc,+G||)2a2]2 

G|| " 



me 



2T:a? ^ 

^ Cikx). 

For the x-direction, the lattice is so infinitely periodic 
that TO, runs from — c» to +c», which should be distin- 
guished from the z-direction. In an actual calculation of 
^{kx) and C{kx), it will be enough to have the first few 
terms of the infinite series because we have a well-defined 
small expansion parameter e~'^ii/". For a further calcu- 
lation of Eq.(pl|), wc feel like getting a simpler form of 
X"^{kz) than in Eq.([10|). If we think of a small damping 
case (kz ^ F), we can simplify Eq.(|lO|) into 

X"'{kz) ~ ;^(fc^)e-i™^-(=^™+=^/2)e''=-('=^"+'=^/2), (15) 

and ImA:^ « r/2|fc2|. On the other hand, we can readily 
perform the remaining integrals for k'^ and fc", 



(18) 

This looks still formidable complicated summations, but 
we can do a systematic summation thanks to a small pa- 
rameter e~'^^l°- (or e~'^\\l°', note cj^ ^ 3c|[). In order to 
retain up to the lowest order of e"'^^' ", we can do a sum- 
mation keeping |to — to'| < 1 and \ml — m"\ < 1. This 
is in the end same order of approximation in S,{kx) and 
Cikx) up to e^'^'^ii/". That is, we have from Eqs.(|l3|) and 

(14) e(fcx)«^[l + 2cos(fc,C||)(c||/a + l)e-^ii/" 



2tt 



+2 cos{2kxC\\ ) (2c|| /a + l)e-2^ii/'' 
+2 co8{3kxC« ) (3c|| /a + l)e-^'=ii/''] 



(19) 



Cik.) 



ICn 



[sin(fc^C|| )e-^ll /" -f 2 sm{2kxC\\)e-'^''« /° 
-t-3sin(3fc^c,|)e~3'=ii/'']. (20) 



We then have three kinds of summations in Eq.(g_8|); we 
define those as a{kz), f3{kz), and 7(fcz)- 



aikz) X (i„e-^^'"'='-'=^) 



/ 



= ^(m — m), (16) 



\ ^ \ ^ \ "^ ~lmkzC±{m-{-m ) ikzC±{m—m ) 
rn m' m" 

x^(to - m')^{m' - m")(l - e^^imfc^c^-j 
1 + (c_L/a + l)2e^2c^/''(l + e-2i'"fe^'=^) 
-|-4cos(fczCj_)(c_L/a -t- l)e 



— c± I a —ln\kz 



+ 2 C0s{2kzCj_){cj_/a + l)2g-2ci/ag-2ImLci ^ (21) 



and also in the same way, 



dkzkzix-ikzTx- (kz) - -, j dkzkz^^j^^^^ 

*''^(TO-TO')e-"^l™-"'l/", 



I3{kz) X (1 - e^2Imfe.c^ 



a a 
Ciin — TO,'). 



\ ^ \ ^ \ ^ —lYnkzC±{m-\-m") ikzC±(m—m") 
m m' m" 

y.C,{m — m)C,(rn ~m )(1 



^ — 21nikzC± > 



(17) 



1 



{c^/afe-^^^^^-{l + e 



— 2lrakzC± ^ 



From Eqs.(Q3)-(17), we can reexpress M(k, i?, w) into a 
form including only the summations for to,'s, 



-iTC0s{2kzCj_){cj_/a) e 



-2c± / a —21mkzC± 



(22) 



-i{k^) X (1 - e" 



-21nikzC± > 



i^^EEE 



— lm.kzC± {m-\-m ) ikzC± (m~m ) 



m m' m" 



x£_{m — Tn)(^{m ^m )(1 — e 



-21mkzC± •) 



■ sm{2k^c±){c±/a + l)(c_L/a)e-^'=^/"e-^^'"'^^'=^ 



sin(fc,cj^)(c_L/a)e-"-^/'^e^^'"'=^"-^. (23) 



Now, all the summation parts of Eq. ( p^ ) are found to be 
real as it should be. Therefore, the photoemission inten- 
sity iCkjEjUi) will be proportional to [— ilmG]^(fcy)] . 
Another factor related to the chain degree of freedom is 
\G2{ky)W= [ReG2(fcj,)]2 + [ImG2(fcy)]2), which gives the 
structure in the final states. Since E + uj ^ E, we can 
assume the final states form a structureless smooth con- 
tinuum. We do not have to make the problem difficult 
unnecessarily. 

Then we can find an analytic solution of the problem. 
The photoemission intensity /(k, E, lu) is 



FIG. 1. Photoemission spectra in the chain direction (|| ky) 
for E II xz. It is seen that the strong fca;-dependencies are in- 
herent in the spectra. In (a), the spectra along (0, 0) ^ (0, ky) 
and in (b), the spectra along (1,0) -^ {^,ky) are provided. 



For an actual calculation and its comparison with ex- 
periment, the reference experiment should be Mizokawa 
et a/.'s pO| . In the experiment, the Fermi energy (the 
width of occupied band) is estimated about 0.2 eV and 
the Fermi wave vector kf = 0.25—. In theory, the 
Fermi energy e/ is given by the linear dispersion within 
Luttinger model such that e/ — vjkj. In the spectral 
function of Luttinger model, i.e., in [— -lmG^(fcy)] , we 
see two features corresponding to the collective charge 
(holon) and spin (spinon) excitations. 



/(k,£:,w) ex 



-ImG+(fcy 



\x{K)\'\xiM'e- 



-lrakzC_i 



X [elC{k,fa{k,) + elkl^(k^fa(k,) 
+el^{k^fl3{k,) + 2e^eyky^{k^)Cik^)aik,) 
+2eyezky^{kxfj{k^) + 2£z£^^(k^)C{k,j:)^{kz)\, 

(24) 



where we note fcz — \ 2[E + uj) — k"^ — k'l 



(a) 

(0.00,0.00) — 
(0.00,0.01) 
(0.00,0 05) 
(0.00,0 10) j 


\ 




- 


, 


^""r^^-i, , 1 


(b) 

(1.00,0.00) — 
(1.00,0 01) 
(1.00,0 05) 
(1.00,0 10) j 

y 


1 

! 

i 

\ 
\ 




- 



-0.9 -0.8 -0.7 -0.6 -0.5 -0.4 



ENERGY (eV) 









- 


(a) 






- 


(1.00,0.00) — 








(1.00,0.05) 






(1.00,0.10) 


\ I 




_^ 3^:i^i^-:- 


— ^— -»=^--= - - -——^ -: 




-0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 

ENERGY (eV) 



FIG. 2. Photoemission spectra for E || 
has a component parallel to the chain. In 



yz, where E 
(a), the spec- 



tra along (1,0) — > i^-yky) and in (b), the spectra along 
(1,—ky) —> {l,ky) are provided. The spectra at (0,0) are 
also given in (a) as a thick solid line. 



The character of features generally rely on the interac- 
tion property through the singularity index a. If a < 1/2, 
there happens a diverging sharp peak at the spinon on- 
set, but if a > 1/2, a converging threshold, a ^ 0.65 is 
estimated from optical study of PrBa2Cu408 ||l^. Here 
we employ this value for PBCO. The charge contribu- 
tion is peaked at —v'f{kf — ky) and the spin excitation at 
—Vf(kf — ky), where we are having the spin-independent 



interaction in mind. 



"/ 



for the holon is taken a little 



larger as v'i/vf ~ 2.5 to be consistent with qualitative 
behaviors of experiment | pO[ . The natural energy unit 
in theory is (Tr/cy)^, whose value is about 5.2 eV from 
C|| ~ 3.8A. From the dispersion, we put ef — Vfkf, 
Vf - 0.15^ and e/ - 0.038(7r/c||)2. The used pho- 
ton energy in the experiment is 29 eV and ~ 145e/, 
from which we take in the calculation the photon en- 
ergy CO ~ 145e/ ~ 5.5(7r/c||)^. However the calculation is 
found quite robust with the parameters. The next thing 
to be considered is a constant damping parameter T. We 
take r — 0.1(7r/c||)^. This corresponds to the mean free 

path A ~ ^ ^ -^^^ ^ 5.3c||. Thinking of c± ~ Scy and 
that the first layer is at ^ , A amounts to governing one 
or two layers. This looks more or less acceptable. In the 
calculation, we take cy/a = 2.5, i.e., c±/a = 7.5. 

Due to the one-dimensionality itself, the chain system 
is expected to show a strong polarization dependency. 
In fact, this is a useful clue in finding one-dimensional 
character submerged in two-dimensional properties. In 
Fig.n^, the photoemission spectra for the photon polariza- 
tion E II xz, more explicitly E — (-^, 0, 4=), are shown. 
In the figure, we attribute the high binding energy struc- 
ture, the sharp peak dispersing between ^ —0.5 eV and 
~ —0.3 eV to the charge excitation (holon) and the low 
binding energy structure, the threshold at ~ —0.2 eV to 
the spin excitation (spinon). When k^ = (Fig.^a)), 
the spinon part is highly suppressed, while when k^ — 1 
(Fig.n^(b)), the spinon part is much enhanced. In the 
experiment pO| , for fc^; = or small k^'s, the spinon 
branch has not been observed and for fc^, ~ 1 , the spinon 
has been so enhanced that it can be observed. Therefore, 
these theoretical and experimental observations has been 
found in good agreement and consistent with each other. 
It need to be noticed that PBCO is an insulator due to a 
charge density wave (CDW) gap opening at ky ^ 0.25 — , 

while the Luttingcr model signifies in principle a metallic 
system. In the study, we provide the results for ky^s away 
from kf because the spin-charge separation is seen only 
for small | ky \ and disappears for ky '^ kf. Further agree- 
ment can be found for E || yz, actually E = (0, -7=, -j=) 
in FigB. In the experiment jlffl, the spectra have been 
suppressed for small kx, but comparable to the configu- 
ration of E II xz when k^ ^ I. This is well understood 
from our theory as shown in Fig.y(a). Comparison of 
the simplest cases (1,0) and (0,0) for both polarizations 
(Figs.|l| and g) shows the fc^, (perpendicular momentum)- 
dependency clearly, whose behaviors are determined by 
the third term in Eq.(|2J), ^{kx)'^l3{kz). It is shown in 
Fig.0 that P(kz) could reinforce or suppress the spinon 
parts depending on the momentum perpendicular to the 
chain, i.e., k^- In the configuration of E || yz, the E 
has a component parallel to the chain, so the symme- 
try of the system, when ky ^> —ky, breaks. Indeed, the 
phase space showing the stronger spinon contribution is 
along (1,0) -^ {1,—ky) rather than (1,0) -^ i^,ky) (for 
ky > 0). Figure ||(b) shows this asymmetry also suc- 



cessfully at least in a qualitative picture. In the cal- 
culations, the asymmetry comes from ky^{kxY^{kz) in 
Eq.(p3). The obtained asymmetry is in fact weaker than 
the experiment and rapidly fades away as |fcj,| increases. 
In the figures, the absolute intensity variations with re- 
spect to several k's would not be very reliable here be- 
cause of many simplifications. Nevertheless, our model 
of the 1/4-filled three-dimensional chain array is found 
to explain very well the one-dimensional behaviors, espe- 
cially the relative charge and spin intensities, reflected in 
the recent ARPES data. 
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FIG. 3. The ratios of a{k^), /3(fcz), and 'yikz) for two sim- 
le cases (1,0) and (0,0) 
and same as for /3 and 7. 



pie cases (1,0) and (0,0) are given, i.e., a = ^[fc^|fc!];^o|fc^^o) 



Another one-dimensional insulator SrCu02 also has 
Cu-0 chains. Even if there is a finite interchain cou- 
pling, it is an order of magnitude weaker than the intra- 
chain coupling, which makes SrCu02 a one-dimensional 
compound. Aside from several details, we find that the 
present calculations should be relevant to the SrCu02 
system as well. In fact, the ARPES experiment on 
SrCu02 shows that the qualitative behaviors are very 
similar to the PBCO case (See Figs. 4 and 5 in Ref. [|3|) 
in that the spinon bands are more apparent in the case 
of large (~ 1) momentum perpendicular to the chain. 



IV. PHOTOEMISSION IN THE SYSTEM OF 
CHAINS : FOR Laa-x-y Ndy SrxCu04 

In the so called stripe phase, there happens the 
one-dimensional spin-charge modulations in the two- 
dimensional Cu02 planes. The stripe phase has attracted 
much attention in connection with its implication to high 
Tc superconductivity since it was proposed from neutron 
scattering in Nd-substituted La2-x-yNdySrxCu04 (Nd- 
LSCO) |12|. It has been also reported that the charge 
transport in Nd-LSCO is observed to be one-dimensional, 
also consistent with the static stripe phase formation [£1| . 
Similar stripe signatures are observed in La2-xSrxCu04 
(LSCO) p3| despite of their dynamical natures, while 
the stripe phase in Nd-LSCO is considered static. In 
this section, we extend the model described in the pre- 



vious section and apply to the Nd-LSCO system. It is 
one of the key issues whether the stripe phase is intrinsi- 
cally metallic or insulating. Here, however, our study will 
proceed based on the simple assumption of metallic one- 
dimensional system. As assumed before, our model has 
completely suppressed the interchain (interstripe) inter- 
action and out-of-chain excitation and may be more suit- 
able for describing the static stripe phase. So Nd-LSCO 
will be preferred in our theory. Nd-LSCO has character- 
istic stripe (one-dimensional) phases, where the stripes 
in adjacent planes are rotated by ^, whose schematic 
sketch is given in Fig.^. It should be noted, however, 
that the ARPES study of LSCO with dynamical stripes 
shows similar characteristics such as one-dimension-like 
low energy electronic structure and the suppression of 
nodal states [ESl. Therefore, the studies of static stripes 
would help to understand the dynamical stripes, too. 




FIG. 4. Schematic view of the static stripe phase observed 
in Nd-LSCO. In the study, the stripe would be modeled as 
a linear charged chain. In an actual experimental situation, 
both configurations have been equally mixed. 
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where i — 1,2 denote the Green function at E and E + oj, 
respectively. The surface is assumed to be at z = 0. The 
Green's function matrix element (klGilk') is 



{k\G,\k') ^ G,{ky)S{ky - k'y) 



X X(fcx)x*(fcx) 



■^<5(fc,-fc; + G||) 



E 

m=0 



x"^ 



G,{k^)5{k, - k'^) 



X X{ky)x'{k'y) 



-i^^(fc,-fc; + G||) 



(k|G,i|k')-l-(k|G,2|k'), 



(26) 



where Gii = 2n7r/c|| is the reciprocal lattice vector for an 
in-plane direction. Following the same line of approxi- 
mation as in the previous PBCO case, we note 



r^ik.) 



Jk,{2c^m+c^/2) 



In the study, we are bearing Lai.28Ndo.6Sro.i2Cu04 in 
our mind as being the reference compound, available in 
the recent ARPES experiment [Q. In this doping, each 
charged stripe is separated by three undoped stripes, so 
the system has a periodicity of Al (l is the lattice constant 
in the in-plane direction), whereas the stripe itself is 1/4 
filled. Then the in-plane interchain (between charged 
stripes) distance c^\ would be cy = 41 and the out-of- 
plane distance would be c_l ~ 2L Here, the wave vector 
unit should be ir/l. In the unit, the Fermi wave vector 
is 1/4. In principle, the starting point is exactly same 
as in PBCO, in the previous section HI except that each 
stripe is assumed a Fermi liquid rather than a Tomonaga- 
Luttinger liquid. For describing the stripes perpendicu- 
larly crossed layer by layer, we write the Green's function 
of the system 



G,;(r,r') = G,(2/-y') ^ x{x - c\im)x* {x' - c^^m) (25) 

rn— — oo 

oo 

X ^ X{z- 2cj_m' - -^)x*{z' - 2cj_m' - -y) 



m'=0 



2 /-^ V ^ 2 

G,{x~x') ^ x(2/-C||m)x*(y'-C||m) 



x"(fc.) 



2/a 



ik^{2c^m+cj_/2) 



(27) 



Then we see the photoemission intensity /(k, E, ui) is ex- 
pressed as (adopting the simplified notation), for A = 



Iik,E,Lo)^G+^AG+,AG^, 

+ G22A.G-^iAG22 

+ G21AGIIAG22 
+ G21AG12AG22 



G'22AG22AG22 

G21AG12AG21 
G22AGj^AG2i 

G22AGj2^^21' 



(28) 
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where we readily recognize only the first two terms 
{G^-^AG^j^AG2i and G^2^Gi2^'^22) are the direct term 
for each y— and x— directed chain set with the perpendic- 
ular interchain distance 2c±, while the others are given 
by the complex interferences of two perpendicular sets. 
The calculation can be done in the same way as in the 
section HI, but its algebra is quite long and tedious. The 
final expression has a slight asymmetry in interchanging 
X— and 2/— coordinates depending on which set is the first 
layer. Thanks to the remaining symmetries in Eq.(^8|), 
nevertheless, it will be enough to calculate the half of 
terms and the other half could be obtained by interchang- 
ing kx <-^ ky and Sx <-> £y if only one account for some 
damping terms. 



In Eg . ( psi) , G21AG11AG21 is exactly same calculation 
as in Eg. (|24|) except for a little changes like cj_ — > 2cj_ 
in a{kz), (3{kz), and 7(^2). Here for a later use, we re- 
define a{kz), I3{kz), "i(kz), S.i'm), and Cip^-) by putting 
c_L -^ 2cj_. For G'j2^Gii^G22 i'^ Eg.(p8|), we have the 
expression, apart from irrelevant constants, 



G+AG+AG 
1 



22 



(29) 
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+4'?(^a)^"(™ - m' + -)e(™' - 771" - -) 
+el^{ky)am -m' + i)C(m' - m" - \) 
+2exeykxC(ky)^(m - m' + -)^(m' - 777" - -) 



1, 



1, 



+eyezC{ky)S,{m - m' + -)C,{m' - m" - -) 
+£v£zC,{ky)C{m -m + -)l{m' - m" - -) 
+£z£xkxi{ky)C,(m - ni + -)^{m' - m" - -) 



+£z£xkxi{ky)£.{m ~ m' + -)C{m' - 777" - -)], 

where r]{k) is another Gii -summation not appeared in sec- 
tion III, 

^(A. + G||f|x(fc + G|()p 



:7^^Ve''=^«"|-C|iN|/a+l|e 

m 
1 C|| 

27ra a 



]^lg-C|||rra|/a 

[l-2cos(fcc||)(c||/a-l)e-=ii/'' 



-2cos(2A:cji)(2cji/a - l)e"2'=ii/"] 
= Vik), (30) 

and the other sophisticated summations for 777,, 777', and 
m" (all r77's > 0) can be managed by considering the 
small parameter e~'^^^°', 



-2IinA;^c_|_ {m-{-7n") i2k^c± {m—m") 



fn m' rn" 



X ^(?71 — ?72 ± — )^(^^ — m =F ~) 
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- sm(2k,ci_)(c±/a)(c±/a + xy-'2c^la^-2\vak,c^ 
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y ^ ^ —4:luikzC±m 



(35) 



The calculated a{kz), P{kz), and "f{kz) are correct to up 
to e'^"^/". Then Gj2^G]^i^G'22 is written down 



^22^^11^*-^22 

--ImG+(fc;, 



ix(fc.)rix(fc.)rix(fc.)re 



2 —SlmkzC 



+el£,{ky)(3{kz) + 2ex£ykxC{ky)a{kz) 
+2eye,C{ky)i{kz) + 2£^£,fc,C(fcj,)7(fc^)]. (36) 

In the calculations of G^j^AGj^j^AGji and G^2^^ii^'-'22' 
we have neglected \G2{ky)\^ (for the latter, \G2{kx)\^) 
based on assuming that the final state distribution may 
be structureless and smooth from E + oj ^ E. Next we 
go to the evaluation of GJ^AGJ^ AG22 + G^^^G^-^AG^i- 
Here we need a more drastic assumption on behav- 
iors of G2{ky), i.e., ImG2(A:j,) ~ 0. Thus we assume 
Gj(/cj,)G^(fca;) or Gj(fca;)G^(fcj,) could be also taken 
just as a positive constant and further |G2(fcjy)p « 
|G2(fc.)|2 « G+{ky)G2{kx) « Gt{kx)G2{ky). Then 
G2+1AG+ AG2-2 is 
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+£y£zkyi(kx)C,{m - m')^{m' - m" - -) 
+£z£xkx£.{kx)({m - ■m')i{m' - m" - -) 
+£z£xC{kx)^{'m - m')l{m' - m" - -)], 
and G22^G^i^G2i is also similarly 
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We can also define a{kz), ${kz), 7i(fcz), and jijkz) as 
the summations for to's relevant in Eqs.(p7[) and (pq); 
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The above a(fcz), f3{kz), 7i(fcz), and 71 (/c^) are also cor- 
rect up to e^^'^^/"'. Now we have done the evaluation of 
G2+1AG+ AG22 + G2+ AG+ AG^l. Using Eqs.(|9|)-(|2l), 
we have 



G2iAG^jAG22 + G22AGj;^AG22 



-ImGf(fcj, 



\x{kxr\x{kyr\x{kz)Ve 



2^-2Imfc^c_L 



X [elkxC{kx)a{kz) + £lklS,{kx)a{kz) 
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+£y£ zkyi{kx)^2{kz) + £z£xkx^{kx)j2{kz) 

+£z£xC{kx)7i{kz)]- 
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(43) 



The other four terms out of Eq. (^8|) can be obtained using 
the symmetry above without actual calculations. From 
Eq. (p8|) , we note the other terms are 



G+ AG+ AG7, = 6-21'"'=-'=^ 



X [G+AG+AG^l; x ^^ y], 



(44) 



where the notation [G^j^AGj^AG^j^; x <;=^ y] means the 
same expression as GjiAGj^AG^j^ with only kx <-> ky, 
£x '^ £y exchanged. When keeping the same notation, 



G2+ AG+ AG^l = e^i-'^^^- 

X [G+AG+AG22; a.^=^y], (45) 

G2iAGj2^^22 + '-'22^^12^^21 

= [G+ AG+ AG22 + G2+ AG+ AG2I ; x<=>y]. (46) 

The total photoemission intensity /(k, E, u) is the sum 
of permutativc contributions (Eq.(28)) obtained above. 

In the experiment of Zhou et a/.'s ||lj], the spectral 
weight distribution has been represented by integrating 
/(k, E, uj) for finite energy windows. From the weight 
distribution, they suggest the superposition of two per- 
pendicular one-dimensional charged stripes for the un- 
derlying electronic structure. Other features such as 
the weight suppression along the d-wave nodal directions 
k = (±1,±1) including the F point does not seem eas- 
ily explained, for which they propose a gap formation 
around \kx\ ~ 1/4 and \ky\ ~ 1/4. The origin of the 
gap, however, does not look very clear. Thus there may 
be left some rooms for another possibility. In our cal- 
culation, the integration of /(k, E, uj) with respect to E 
can be readily done. From Eqs.(^), (p5|), and (^), the 
S-dependency is in [— ilmGf (fcj,)] or [— j^lmGf (k^)] 

< 



Ky, 



\E\ 



and through kz- For k^ — \ r2.{E + uj) — k 

0.1 eV from the band width, while w ~ 0(10) eV in 
actual experiments. So we could approximate kz ~ 



2uj~kl 



ky. Then the remaining integral is special, 



dE 



--lmG+{ky(x))I{E) 



nk„ 



(47) 



where f{E) is the Fermi-distribution function. 
\—-IniGj^{ky(^x))] is the spectral function of the stripe, 
which determines its band structure. Because our start- 
ing point is the noninteracting chain (stripe), we should 
assume a suitable band dispersion relevant to the one- 
dimensional stripe. The corresponding band by the 
model calculation pj] from vertical stripes is very flat 
especially near (7r,0), which is believed to have the one- 
dimensional character in the complicated real band dis- 
persion. In the experiments |25|, its band width is es- 
timated about 30 meV and should be the integration 
energy window in our study. So we assume the band dis- 
persion of the stripe to be bky — bki for y-directed stripe, 
where 6 ~ 0.18 from {Tr/lf ~ 2.64 eV {I ~ 5.34A). It 
should be noticed that, despite one-dimensional prop- 
erties in the Nd-LSCO system, there have been no 
compelling evidences for Tomonaga-Luttinger-like elec- 
tronic structure (i.e., such as the spin-charge separation) 



in the stripe phase. Therefore, even if the matrix ef- 
fects are explored based on one-dimensional stripes, we 
assume the Fermi-like momentum distribution nt , , , 
which is motivated by our interests in the polariza- 
tion, photon energy and momentum dependencies of the 
ARPES intensity, not in the energy distribution curves 
(EDO's). Then we get the expression /(k, cj), defining 
I{k,uj) = JdEI{k,E,u;)f{E), 

/(k,c.)«nfeJx(fc.)nx(fc.)Pe-i-^^^- 

X [slak.fc^ikz) + elklak.)^c.{kz) 

+el^{kxff3{kz) + 2exeykyS,{kxK{kx)a{kz) 
+2eyezky^{kxfj(kz) + 2ezex^{kx)C{kxh{kz)] 



-2Imfe^c 



X [kx 



f^y 1 ^x 



+ 2n,Jx{kx)\'\x{ky)\'\x{kz)\'e-''^~''^- 
X [elkxC{kx)a{kz) + elkl^{kx)a{kz) 

+el£_{kx)l3{kz) + exeyky({kx)a{kz) 
+SxSykxky^{kx)a{kz) + Sy£zky^{kx)ji{kz) 
+£y£zkyi{kx)^2{kz) + £z£xkxS.{kx)l2{kz) 

+£z£xC{kx)ll{kz)] 



[r^x 



£«] 



+ nkMk^)?\x{ky)\'\x{kz)?e-^'^^'^-^- 
X [£lkU{ky)a{kz) + elT]{ky)a{kz) 
+el£_{ky)P{kz) + 2exeykxC{ky)a{kz) 



2ey£zC{ky)l{kz) + 2£zexkxS,{ky)^{kz)] 

>ky,£x^^=> £y\, (48) 



g2Imfc.c^ X [kx 



where [kx <J=^ ky,£x <J=> Sy] denotes the term right be- 
fore it with kx *-> ky and Sx <-> £y exchanged. According 
to Eq. ( |4q ) , because of F (or Im^^ ) , we expect to get an 
asymmetry in the weight distribution (with respect to kx 
and ky) depending on which is the first layer because the 
mean free path could govern only the first few layers in 
realistic situations. To compare with experiment jlj], we 
average the spectra with respect to the stripe direction, 
i.e. for the two configurations in Fig.0. 

From Eq.(48), we evaluate the spectral weight distri- 
bution for E ± xy and provide the results in Fig||. Note 
the momentum distribution Uk f^-. readily gives two sets 
of Fermi surfaces defined by [kx\ = 1/4 and [ky[ — 1/4. 
In the figure, we have found two kinds of distributions de- 
pending on the photon energy. Figures o(c) and H(d) may 
be understood just as the superposition of two perpen- 
dicular sets, whereas Figs.^(a) and||(b) are very interest- 
ing. As in the experiment, we can observe the spectral 
weight reduction along the d-wave nodal direction and 
also around the F point, and simultaneously the spectral 
confinement around (0, ±1) and (±1,0). This must be 
the geometric (interference) effect raised by the dipole 
matrix in our study. More clearly, when Ex = Sy = 0, 
/(k, E) is determined by the functions P{kz), f3{kz), and 
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f3{kz), of which f3{kz) must be dominant. Thus we can 
see that /3(fc^) is mainly attributed to !{]<., E), that is, 
G22^GtiAG22 + G2iAG^2^G2i, which are the interfer- 
ence terms, not direct terms. If we try to explain the 
experiment based on the present results, a gap opening 
is not required and metallic stripes look enough to repro- 
duce the observed weight distribution. Unfortunately, 
however, it is also found in our study that two types of 
spectral distributions appears alternatively depending on 
the photon energy. Experimentally, the intensity confine- 
ment near (0, ±1) and (±1, 0) is robust when uj ^ 20 eV, 
30 eV, and 50 eV. In FigJ6| we give the spectral intensities 
at the r point for the photon energy cu. 
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FIG. 5. Spectral weight distributions for the linearly polar- 
ized photon, especially perpendicular ((a) and (c)) or nearly 
perpendicular ((b) and (d)) to the layer are given. We have 
two kinds of distribution depending on the photon energy. 
(a) and (b) are for to — 15(7r//)^ and (c) and (d) are for 
CO = 20(7r//)^. The arrow in the figures ((b) and (d)) denotes 
the direction of the in-plane electric field, l/a = 2.5 is used. 
The wave vector unit is n/l. 



In Fig.0, the spectral distribution for E || xy is given, 
which is quite a different distribution from Fig.ra. In the 
figure, the intensity at the F point is suppressed and there 
happen to be strong weights on the diagonals, i.e., (0, 0) 
to (±1,±1). According to the recent experiment p5[ , 
for E parallel to the surface, weak intensities appear also 
in the diagonal directions in addition to the intensities 
at (0,±1) and (±1,0), which would not be observed or 
observed weaker in the case of E perpendicular to the 
surface. These additional complex structures have been 
suggested as the overlap of two-dimensional characters 
from stripe fluctuations on the simple one-dimensional 



stripe picture, namely as part of the two dimensional 
Fermi surface. On the other hand, such fluctuations were 
not considered in our study. 




PHOTON ENERGY 



FIG. 6. Spectral intensities at the F point are given with 
respect to the photon energy when E _L xy. The energy unit 
is in/l)\ 
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FIG. 7. The spectral weight distribution for E 
— 15(7r//)^. The arrow is the direction of E. 
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As shown in the figures, we do not stress the re- 
stricted agreement much. The stripe system, Nd-LSCO 
is a more complicated system compared to the chain sys- 
tem, PBCO. The real electronic structure may consist of 
both one-dimensional and two-dimensional components. 
This is very important and can explain why the present 
approach is good for PBCO, while less satisfactory for 
Nd-LSCO. The photon energy behavior in Figsjg and ||, 
however, provides outlook for further studies along the 
present line. The first thing may be to go beyond the 
primary plane wave assumption for the photoelectron. 
The photoelectron wave function should be in principle 
calculated using the scattering theory including the sur- 
face. Also the excitation perpendicular to the stripes, 
neglected here, would introduce the surface Green's func- 
tion Q. Wc would then hopefully expect more correct 
transition matrix effects, and also more robust photon en- 
ergy and polarization dependencies of the spectra. How- 
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ever, we do not mean here the matrix effect would be 
enough to understand experiments. We have regarded 
stripes as a cohection of noninteracting metaUic chains, 
but the interaction among the chains as weU as the inter- 
action between the chain and the two-dimensional back- 
ground electronic structures may play a role essentially 



V. CONCLUSIONS 

We have studied the dipole matrix effects of the pho- 
toemission in the system of linear metallic chains for 
two characteristic configurations. Assuming the chains 
(stripes) are perfectly one-dimensional by neglecting the 
transverse excitations and interchain interactions, we 
could write down the Green's function of the total sys- 
tem. Pendrey's formula enables us to manage the dipole 
matrix effects from knowledges of the single-particle 
Green's function. The calculations are very well con- 
trolled by small parameters e~'^\\l°- or e""^^/". 

Interesting k- and E-dependencies in the photoemis- 
sion spectra in PBGC system have been observed, which 
would be thought of as the transition matrix effects. We 
have evaluated the spectral intensities accounting for the 
dipole matrix in the system of chains based on the Lut- 
tinger model. PBGO is found to have one-dimensional 
chains which behave like Tomonaga-Luttinger liquid by 
recent ARPES studies. For the photoemission formula, 
we have adopted the Pendrey's formula. The relative in- 
tensities of separated liolons and spinons depending on 
the wave vector k (anisotropics), or the photon polar- 
ization (linear polarization, E) obtained in the study are 
found consistent with the recent experiments. It is also 
found that the present study would be applicable to an- 
other system having insulating Gu-0 chains, SrGu02. 

We have also considered the system more extended 
than PBGO. We have had interests in the static stripe 
phases observed in Nd-LSGO system. In experiments, 
their one-dimensional characters are observed from the 
spectral distribution in k-space, but not perfectly under- 
stood from the stripe picture alone. We have studied its 
deviations arising from the dipole matrix effects. In the 
same way as for PBGO, we have calculated the integrated 
spectral distribution. We then found a distribution quite 
consistent with experiment (for E _L xy), where the in- 
tensity around the P point or along the d-wave nodal 
directions is strongly suppressed. We have also found, 
however, another distribution showing a strong accumu- 
lation of intensities around the P point depending on 
the photon energy. In the result, two types of distri- 
butions appear alternatively as a function of photon en- 
ergy. Such more or less shaky behaviors may be due to 
the simple plane wave assumption for the photoelectron 
state. For E || xy^ additional complex structures are 
more apparently observed in experiment. It is insisted 
that the system should consist of two components, both 



one-dimensional and two-dimensional characters even in 
the static stripe phase in Nd-LSGO [^. E-dependent 
behaviors in our calculations, however, come from the 
interferences of noninteracting stripes (chains) and may 
be expected not fully justified in the stripe system. In the 
stripe system much two-dimensional characters induced 
by stripe fluctuations and interactions of the stripes with 
the two-dimensional background exist unlike a pure chain 
system. This is the reason that our study is successful 
for PBGO, whereas less successful in Nd-LSGO. 

Finally, we may need to clarify the scope or limit for 
the present photoemission approach. We have thought 
of the noninteracting (noninteracting among chains and 
also between chains and the rest of the system) one- 
dimensional chains (stripes) for both PBGO and Nd- 
LSGO. Therefore, in such a sense, the matrix effects ad- 
dressed in our study are purely geometric. This must 
be quite drastic approximations, but can be well con- 
sistent with our original purpose of studying the matrix 
effects. Further investigations for the dipole matrix effect 
are suggested to take into account better considerations 
for the stripes (allowing the transverse excitations) and 
thus leading to the surface Green's function, which would 
make it possible to describe the dynamical stripe phase 
too. We should also think of how good the plane wave for 
the photoelectron could be in these anisotropic systems. 
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